In situ stress field of eastern Australia     

R. R. Hillis  J. R. Enever  S. D. Reynolds 《Australian Journal of Earth Sciences》2013,60(5):813-825

New in situ data based on hydraulic fracturing and overcoring have been compiled for eastern Australia, increasing from 23 to 110 the number of in situ stress analyses available for the area between and including the Bowen and Sydney Basins. The Bowen Basin displays a consistent north‐northeast maximum horizontal stress (σ_H) orientation over some 500 km. Stress orientations in the Sydney Basin are more variable than in the Bowen Basin, with areas of the Sydney Basin exhibiting north‐northeast, northeast, east‐west and bimodal σ_H orientations. Most new data indicate that the overburden stress (σ_V) is the minimum principal stress in both the Bowen and Sydney Basins. The Sydney Basin is relatively seismically active, whereas the Bowen Basin is relatively aseismic. Despite the fact that in situ stress measurements sample the stress field at shallower depth than the seismogenic zone, there is a correlation between the stress measurements and seismicity in the two areas. Mohr‐Coulomb analysis of the propensity for failure in the Sydney Basin suggests 41% of the new in situ stress data are indicative of failure, as opposed to 13% in the Bowen Basin. The multiple pre‐existing structural grains in the Sydney Basin further emphasise the difference between propensity for failure in the two areas. Previous modelling of intraplate stresses due to plate boundary forces has been less successful at predicting stress orientations in eastern than in western and central Australia. Nonetheless, stress orientation in the Bowen Basin is consistent with that predicted by modelling of stresses due to plate boundary forces. Variable stress orientations in the Sydney Basin suggest that more local sources of stress, such as those associated with the continental margin and with local structure, significantly influence stress orientation. The effect of local sources of stress may be relatively pronounced because stresses due to plate boundary forces result in low horizontal stress anisotropy in the Sydney Basin.  相似文献   

Hydraulic diffusivity and macrodispersivity calculations embedded in a geographic information system     

F.A.L. Pacheco 《水文科学杂志》2013,58(4):930-944

Abstract

A numerical technique is presented whereby aquifer hydraulic diffusivities (D) and macrodispersivities (α) are calculated by linear equations rewritten from flow and solute transport differential equations. The approach requires a GIS to calculate spatial and temporal hydraulic head (h) and solute concentration gradients. The model is tested in Portugal, in a semi-confined aquifer periodically monitored for h and chloride/sulphate concentrations. Average D (0.46 m²/s) and α (1975 m) compare favourably with literature results. The relationship between α and scale (L) is also investigated. In this context, two aquifer groups could be identified: the first group is heterogeneous at the “macroscopic” scale (solute travelled distances <1 km), but homogeneous at the “megascopic” scale. The overall scale dependency in this case is given by an equation of logarithmic type. The second group is heterogeneous at the macroscopic and megascopic scales, with a scale dependency of linear type.

Citation Pacheco, F.A.L., 2013. Hydraulic diffusivity and macrodispersivity calculations embedded in a geographic information system. Hydrological Sciences Journal, 58 (4), 930–944.  相似文献   

Using a H∞ filter assimilation procedure to estimate root zone soil water content     

Haishen Lü  Xiaoli Li  Zhongbo Yu  Robert Horton  Yonghua Zhu  Zhenchun Hao  Long Xiang 《水文研究》2010,24(25):3648-3660

Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Influence of pore size and geometry on peat unsaturated hydraulic conductivity computed from 3D computed tomography image analysis     

F. Rezanezhad  W. L. Quinton  J. S. Price  T. R. Elliot  D. Elrick  K. R. Shook 《水文研究》2010,24(21):2983-2994

In organic soils, hydraulic conductivity is related to the degree of decomposition and soil compression, which reduce the effective pore diameter and consequently restrict water flow. This study investigates how the size distribution and geometry of air‐filled pores control the unsaturated hydraulic conductivity of peat soils using high‐resolution (45 µm) three‐dimensional (3D) X‐ray computed tomography (CT) and digital image processing of four peat sub‐samples from varying depths under a constant soil water pressure head. Pore structure and configuration in peat were found to be irregular, with volume and cross‐sectional area showing fractal behaviour that suggests pores having smaller values of the fractal dimension in deeper, more decomposed peat, have higher tortuosity and lower connectivity, which influences hydraulic conductivity. The image analysis showed that the large reduction of unsaturated hydraulic conductivity with depth is essentially controlled by air‐filled pore hydraulic radius, tortuosity, air‐filled pore density and the fractal dimension due to degree of decomposition and compression of the organic matter. The comparisons between unsaturated hydraulic conductivity computed from the air‐filled pore size and geometric distribution showed satisfactory agreement with direct measurements using the permeameter method. This understanding is important in characterizing peat properties and its heterogeneity for monitoring the progress of complex flow processes at the field scale in peatlands. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

A finite element approach to the simulation of hydraulic fractures with lag     

Michael J. Hunsweck  Yongxing Shen  Adrian J. Lew 《国际地质力学数值与分析法杂志》2013,37(9):993-1015

We presented a finite‐element‐based algorithm to simulate plane‐strain, straight hydraulic fractures in an impermeable elastic medium. The algorithm accounts for the nonlinear coupling between the fluid pressure and the crack opening and separately tracks the evolution of the crack tip and the fluid front. It therefore allows the existence of a fluid lag. The fluid front is advanced explicitly in time, but an implicit strategy is needed for the crack tip to guarantee the satisfaction of Griffith's criterion at each time step. We enforced the coupling between the fluid and the rock by simultaneously solving for the pressure field in the fluid and the crack opening at each time step. We provided verification of our algorithm by performing sample simulations and comparing them with two known similarity solutions. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Numerical modeling of fluid‐driven fractures in cohesive poroelastoplastic continuum     

E. Sarris  P. Papanastasiou 《国际地质力学数值与分析法杂志》2013,37(12):1822-1846

In this article, we investigate the main parameters that influence the propagation of a fluid‐driven fracture in a poroelastoplastic continuum. These parameters include the cohesive zone, the stress anisotropy, and the pore pressure field. The fracture is driven in a permeable porous domain that corresponds to weak formation by pumping of an incompressible viscous fluid at the fracture inlet under plane strain conditions. Rock deformation is modeled with the Mohr–Coulomb yield criterion with associative flow rule. Fluid flow in the fracture is modeled by the lubrication theory. The movement of the pore fluid in the surrounding medium is assumed to obey the Darcy law and is of the same nature as the fracturing fluid. The cohesive zone approach is used as the fracture propagation criterion. The problem is modeled numerically with the finite element method to obtain the solution for the fracture length, the fracture opening, and the propagation pressure as a function of the time and distance from the pumping inlet. It is demonstrated that the plastic yielding that is associated with the rock dilation in an elastoplastic saturated porous continuum is significantly affected by the cohesive zone characteristics, the stress anisotropy, and the pore pressure field. These influences result in larger fracture profiles and propagation pressures due to the larger plastic zones that are developing during the fracture propagation. Furthermore, it is also found that the diffusion process that is a major mechanism in hydraulic fracture operations influences further the obtained results on the fracture dimensions, plastic yielding, and fluid pressures. These findings may explain partially the discrepancies in net pressures between field measurements and conventional model predictions. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

Simulation of hydraulic fracture propagation near a natural fracture using virtual multidimensional internal bonds     

Z. Zhang  A. Ghassemi 《国际地质力学数值与分析法杂志》2011,35(4):480-495

A virtual multidimensional internal bond (VMIB) model developed to simulate the propagation of hydraulic fractures using the finite‐element method is formulated within the framework of the virtual internal bond theory (VIB) that considers a solid as randomized material particles in the micro scale, and derives the macro constitutive relation from the cohesive law between the material particles with an implicit fracture criterion. Hydraulic pressure is applied using a new scheme that enables simulation of hydraulically driven cracks. When the model is applied to study hydraulic fracture propagation in the presence of a natural fracture, the results show the method to be very effective. It shows that although the in situ stress ratio is the dominant factor governing the propagation direction, a natural fault can also strongly influence the hydraulic fracture behavior. This influence is conditioned by the shear stiffness of the fault and the distance to the original hydraulic fracture. The model results show that when the fault is strong in shear, its impact on hydraulic fracture trajectory is weak and the hydraulic fracture will likely penetrate the fault. For a weak fault, however, the fracture tends to be arrested at the natural fault. The distance between the fault and the hydraulic fracture is also important; the fault influence increases with decreasing distance. The VMIB does not require selection of a fracture criterion and remeshing when the fracture propagates. Therefore, it is advantageous for modeling fracture initiation and propagation in naturally fractured rock. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

On the hydraulic intake shape factor for a circular opening located at an impervious boundary: Influence of inclined stratification     

A. P. S. Selvadurai 《国际地质力学数值与分析法杂志》2011,35(6):639-651

In this paper, we develop an exact closed form solution for a circular entry point located at the interface between an impermeable material and a stratified porous medium, where the principal plane of hydraulic isotropy is inclined to the interface. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Soil water movement under a bottomland hedgerow during contrasting meteorological conditions     

R. Ghazavi  Z. Thomas  Y. Hamon  Ph. Merot 《水文研究》2011,25(9):1431-1442

  相似文献   

Plant drought tolerance trait is the key parameter in improving the modeling of terrestrial transpiration in arid and semi-arid regions   总被引：1，自引：0，他引：1       下载免费PDF全文

Xintao Liu  Xingjie Lu  Shulei Zhang  Zhongwang Wei  Nan Wei  Shupeng Zhang  Hua Yuan  Wei Shangguan  Shaofeng Liu  Jianfeng Huang  Lu Li  Xiulan Ye  Jinxuan Zhou  Wenke Hu  Yongjiu Dai 《大气和海洋科学快报》2022,15(1):35-41

陆面过程蒸腾作用的模拟制约着天气,气候降水预测的精确度.近几十年来,为了更好地描述植被蒸腾的水力约束,陆面过程模式发展了基于植物性状的植物水力胁迫方案.然而,我们对于植物性状在蒸腾模拟中的地位仍然缺乏了解,植物性状对蒸腾的重要性仍需进一步量化.本研究利用Morris方法评估植物性状参数在通用陆面模式植物水力胁迫方案(CoLM-P₅₀HS)中的重要性,针对17种植物性状,筛选出最为重要的:耐旱性状(P₅₀),气孔性状,和光合作用性状.在12个FLUXNET站点中,参数的重要性由归一化敏感度来衡量.P₅₀的重要性随着降水的减少而增加,而气孔性状和光合作用性状的重要性则随着降水的减少而减少.在干旱或半干旱地区,P₅₀比气孔性状和光合作用性状更重要,这意味着当植物经常经历干旱时,水力安全策略比植物生长策略更关键.而耐旱性状的巨大变异性进一步暗示了多种植物水力安全策略的共存.忽视P₅₀的变异性可能会对陆面过程模式蒸腾作用的模拟造成严重误差.因此,为了更好地表示植物水力功能的变异性,需要增加对耐旱性状的观测并耦合到陆面模式中.  相似文献